Effects of mirror dark matter on neutron stars 
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If dark matter is made of mirror baryons, they are present in all gravitationally bound structures. 
Here we investigate some effects of mirror dark matter on neutron stars and discuss possible ob- 
servational consequences. The general-relativistic hydrostatic equations are generalized to spherical 
objects with multiple fluids that interact by gravity. We use the minimal parity-symmetric exten- 
sion of the standard model, which implies that the microphysics is the same in the two sectors. We 
find that the mass-radius relation is significantly modified in the presence of a few percent mirror 
baryons. This effect mimics that of other exotica, e.g., quark matter. In contrast to the common 
view that the neutron-star equilibrium sequence is unique, we show that it depends on the relative 
number of mirror baryons to ordinary baryons. It is therefore history dependent. The critical mass 
for core collapse, i.e., the process by which neutron stars are created, is modified in the presence of 
mirror baryons. We calculate the modified Chandrasekhar mass and fit it with a polynomial. A few 
percent mirror baryons is sufficient to lower the critical mass for core collapse by ~ O.I M0. This 
could allow for the formation of extraordinary compact neutron stars with low mass. 

PACS numbers: 11.30.Er, 26.60.Kp, 95.35.+d, 97.60.-s 
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I. INTRODUCTION 

Observations indicate that non-luminous matter con- 
tributes a significant part of the mass in galaxies. This 
"dark matter" (DM) cannot be made of known particles, 
but is of unknown type. In this paper we describe some 
potential effects of DM on neutron stars (NS), which can 
be tested with future observations of NS. These effects 
could be important also for the interpretation of the ob- 
servational results. If DM exists, it should be present in 
all astrophysical objects, unless there is an efficient and 
unexpected segregation mechanism. DM can be present 
already in the formation processes and it can be accreted 
subsequently from the environment. Both of these pro- 
cesses should naturally occur, but the amount of DM 
that is trapped inside different objects should depend on 
the nature of the DM particles, the type of objects and 
their history. In general we expect that critical phases 
of stellar evolution could be modified in the presence of 
relatively small amounts of DM. 

The most studied hypothesis is that DM is made of 
weakly interacting massive particles (WIMPs), which 
naturally appears in supersymmetric extensions of the 
standard model (SM) of particle physics. Should they 
exist, WIMPs accumulate in NS due to elastic scatter- 
ing on nucleons. If the mass of accreted WIMPs reaches 
a critical value they become self-gravitating and form a 
dense core supported by degeneracy pressure. The dark 
core collapses into a black hole if it reaches the Chan- 
drasekhar mass, M ~ Mpi/m'jjj^.j. This limits thepossi- 
ble density, cross section and mass of WIMPs See 
also [3j . The expected high mass of WIMPs gives a low 
Chandrasekhar mass for the DM core, they therefore do 
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not affect the mass and radius of NS significantly. The 
annihilation of WIMPs in the DM core produces heat and 
old NS should therefore have a steady-state temperature 
of about lO"*^ K [3]. The coldest observed NS have tem- 
peratures in the range 10^ — 10^ K. It is observationally 
difficult to detect NS with significantly lower tempera- 
tures, because the luminosity scales as T"^. The effects 
of WIMPs on NS therefore seems to be of little practi- 
cal concern. WIMPs are, however, only one class of DM 
particles and some observations suggest that the model 
is oversimplified, see, e.g., Q and references therein. 

Mirror matter is another viable DM candidate that 
emerges if one, instead of (or in addition to) assuming 
a symmetry between bosons and fermions - supersym- 
metry, assumes that nature is parity symmetric. It is 
a matter of fact that the weak nuclear force is not par- 
ity symmetric. The main theoretical motivation for the 
mirror-matter hypothesis is that it constitutes the sim- 
plest way to restore parity symmetry in the physical laws 
of nature. When Lee and Yang proposed the non-parity 
of weak interactions in 1956, they mentioned also the 
possibility to restoreparity by doubling the number of 
particles in the SM jg]. Thereby the Universe is divided 
into two sectors with opposite handedness that interact 
mainly by gravity. In the minimal parity-symmetric ex- 
tension of the SM 0, Q , the group structure is G (ii G, 
where G is the gauge group of the SM. In this model the 
two sectors are described by the same lagrangians, but 
where ordinary particles have left-handed interactions, 
mirror particles have right-handed interactions. Except 
for gravity, mirror matter could interact with ordinary 
matter via so-called kinetic mixing of gauge bosons, or 
via unknown fields that carry both ordinary and mirror 
charges. If such interactions exist they must be weak 
and we therefore neglect them here. Since photons do 
not interact with mirror baryons, or interact only via the 
weak kinetic mixing, mirror matter constitute a natu- 
ral candidate for the DM in the Universe. The study of 



2 



the cosmological implications of mirror matter is well- 
defined, because the microphysics of mirror matter is the 
same as that of ordinary matter. Many consequences of 
mirror matter for particle physics and astrophysics have 
been studied during the last decades. The reader can re- 
fer to f9] for a review of the history of mirror matter and 
a list of most relevant papers published before 2006. 

Like their ordinary counterparts, mirror baryons can 
form atoms, molecules and astrophysical objects such 
as planets, stars and globular clusters. However, even 
though the microphysics is the same in both sectors, the 
chemical content of the mirror sector should be different, 
because the cosmology must be different. In particu- 
lar. Big Bang nucleosynthesis (BBN) requires that the 
mirror sector has lower temperature than the ordinary 
sector [l^, [HI. This has implications for the thermo- 
dynamics of the early Universe [l^, [l^l and for the key 
cosmological epochs. Analytical results and numerical 
calculations of BBN (Tol . [Tsj show that the mirror sector 
should be helium dominated, and the abundance of heavy 
elements is expected to be higher than in the ordinary 
sector. Invisible stars made of mirror baryons are candi- 
dates for Massive Astrophysical Compact Halo Objects 
(MACHOs)^^hich have been observed via microlensing 
events [13, [iH [iBl- Mirror stars contain more helium 
and less hydrogen, and therefore have different proper- 
ties than ordinary stars [iTj . The consequences of mirror 
matter on primordial structure formation, the cosmic mi- 
crowave background and the large-scale-structure distri- 
bution of matter have been studied [H, [H, [M [lH, [H, [H . 
These studies provide stringent bounds on the mirror sec- 
tor and prove that it is a viable candidate for DM. In 
addition, mirror DM is one of the few potential explana- 
tions for the recent DAM A/LIBRA annual modulation 
signal [23,[2i,[2^. 

The properties of NS are intimately connected to the 
equation of state (EoS) of matter at densities well beyond 
nuclear saturation density, no ~ 0.16 fm~'^. NS there- 
fore are natural laboratories for the exploration of bary- 
onic matter under extreme conditions, complementary 
to those created in terrestrial experiments with atomic 
nuclei and heavy-ion collisions. A stiff EoS at high den- 
sity is needed to explain NS with high mass [27| and 
large Raa 128|]. Heavy- ion collision data for kaon pro- 
duction and elliptic flow provide an upper limit on the 
stiffness of the EoS [1^, [s^]- By combining these con- 
straints it is possible to test and exclude certain models 
of high-density EoS [Slj, [s^ • In this context it is impor- 
tant to be aware of the potential effects of DM on the 
properties of NS. A long-debated question in compact- 
star physics is whether quark matter exist in the core of 
NS and whether there are unambiguous observables as- 
sociated with that. Other exotic states of matter, e.g., 
meson condensates and hyperons could also exist in the 
core. For recent reviews, see [3l,[33,[3l]. For an example 
model of hybrid stars (stars with a quark-matter core en- 
closed in a nuclear-matter shell) that are consistent with 
present observational and experimental constraints, see 



|32l |. A typical consequence of exotica is that the maxi- 
mum NS mass and the radii of NS with typical masses, 
M ^ 1.35 Mq, becomes lower. Constraints on the mass 
and mass-radius relation of compact stars are therefore 
used as indicators for the presence/absence of exotica 
in compact stars. As an example, see the claim in [s^ 
and the counter-examples provided in (37j . In this pa- 
per we show that if mirror matter (or, in principle, some 
other form of stable self-interacting DM) accumulates in 
NS, the equilibrium sequence could be significantly mod- 
ified and the effect is similar to that of traditional exotic 
phases of matter. The NS equilibrium sequence is di- 
rectly related to the ground-state equation of state and 
it is therefore commonly assumed to be unique. We show 
that if DM accumulates in NS this is not necessarily the 
case. See also [ssj . where potential effects of mirror mat- 
ter on NS are discussed qualitatively. 



II. COMPACT STAR SEQUENCE 

The structure of NS with a mirror-matter part is deter- 
mined by hydrostatic equations, which are similar to the 
well-known Oppenheimer-Volkoff (OV) equations. Here 
we repeat some essential steps in the derivation of the OV 
equations, and then we generalize the result to include 
mirror matter. The starting point is the line element, 
dr^ = g^^dx^^dx" , of static isotropic spacetime 



dr' 



and the energy-momentum tensor of a perfect fluid 



-pg^" + {p + p)u^u\ 



(2) 



Here, p is the pressure and p the energy density, which 
includes rest-mass energy. The Einstein field equation 
simplifies to (units are chosen such that G — c — 1) 
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where 



M(r) = 47r / dff'^p{f). 
Jo 



(3) 
(4) 
(5) 

(6) 



The metric should match the Schwarzschild solution at 
the surface of the star, which is located at r = i?. It then 
follows from ([3]) that the gravitational mass of the star 
is M = M{R). Similarly, the differential equation ^ for 
h'{r) is subject to the Schwarzschild boundary condition, 
so the integration constant is fixed. In combination with 
an EoS, p = p{p), equation ^ determines how the pres- 
sure and density varies with r inside a star. Given an EoS 
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and a central pressure, p{r = 0), ([I])-® are integrated 
until the surface is reached, where p{r = R) = 0. By 
varying the central pressure, a one-parameter sequence 
of equilibrium solutions with different M and R is ob- 
tained. 

Newtonian physics applies to systems where rest-mass 
energy dominates. In general relativity, all forms of en- 
ergy are sources of gravity. In particular, the curvature 
of spacetime inside a compact star depends not only on 
the energy-density distribution of matter, but also on the 
pressure, see In the presence of mirror matter the 
metric is affected by the energy content in both sectors 



Pir) = Poir) +PM{r), 
p{r) = po{r) + PM{r). 



(7) 

(8) 



Here, O (M) refers to the ordinary (mirror) sector. The 
metric functions, A(r) and v{r), applies to both sectors 
and are given by ©-(HI) with the replacements ©-([HI). 
Equation ([5]) separates for the two sectors 

^ = -bMM + mM]^, (10) 

because here we assume that particles and mirror parti- 
cles interact by gravity only 51]. This result can be gen- 
eralized to any number of fluids that interact by gravity. 
In principle the two sectors could interact by other forces 
than gravity, e.g.^ by photon-mirror-photon kinetic cou- 
pling, which gives charged mirror particles an effective 
electric charge in the ordinary sector that is about eight 
or nine orders of magnitude lower than that of ordinary 
charged particles [33|- The effect of such eventual weak 
interactions on the equilibrium structure of a compact 
star should be small and we therefore consider only the 
gravitational interaction here. 

Equations ^ and ([6l)- (fT0| are the general-relativistic 
hydrostatic equations for spherical stars with a mirror 
matter part. By choosing different central pressures in 
the two sectors, a star with different radius and baryon 
number in the two sectors is obtained. In Fig. [1] we plot 
compact star sequences with different number of baryons 
in the mirror and ordinary sectors. The EoS used here 
is based on a relativistic mean-field model of nuclear 
matter 'i^, which is combined with the Baym-Pethick- 
Sutherland (BPS) EoS for the crust ^l,]. In this work we 
approximate the EoS of mirror nuclear matter with that 
of ordinary nuclear matter, i.e., we use the same EoS in 
the two sectors. We motivate this approximation in the 
following way: NS are formed from the iron core of stars, 
i.e., of nuclei with maximum binding energy formed at 
the end of nuclear burning. They quickly become cold 
on the nuclear energy scale, because of neutrino emis- 
sion, and matter thereafter is in the ground state. The 
situation is different for the mirror-matter part of an NS. 
Different types of mirror nuclei can be accreted, e.g., from 
a binary companion mirror star, and an eventual mirror 
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FIG. 1: (Color online) Neutron-star sequences for differ- 
ent number of baryons in the mirror sector, Nm/{No + 
Nm) = 0%, 1%, 10%, 25% and 50%. Here, M is the total 
gravitational mass of the star and Ro (Rm) is the coordinate 
radius of the surface in the ordinary (mirror) sector. Here 
we consider 7Vi\/ < No only, but the results can be gener- 
alized to the opposite case by the replacement Ro ^ Rm. 
These sequences correspond to a nuclear relativistic mean- 
field equation of state [4^ and the crust is modeled with 
the BPS equation of state Included in the figure are 

also lines of constant surface redshift (thin dashed lines), 
z = 1/^1- 2M/R - 1, and 
ted lines), Rao — R/ 



"radiation radius" (thin dot- 
2M / R. The latter quantity can 
be constrained via the surface emission from isolated neutron 
stars and is, in addition to the mass, an important indirect 
observable [2^. The surface redshift has not yet been reli- 
ably observed, but with future high-resolution observatories 
and a better understanding of neutron-star atmospheres this 
can become an important observable. 



core of the progenitor star is not necessarily made of iron. 
The mirror core is, however, compressed into a compact 
object by the strong gravitational field of the NS. Mirror 
nuclei are therefore disintegrated into mirror neutrons, 
mirror protons and light mirror nuclear clusters, which 
should equilibrate with respect to weak nuclear reactions, 
just like neutrons and protons do in ordinary NS. In the 
minimal parity-symmetric extension of the SM (see the 
Introduction) the lagrangian of the mirror sector is iden- 
tical to that of the ordinary sector. In this context it 
therefore seems reasonable to approximate the EoS of 
high-density mirror nuclear matter with that of ordinary 
nuclear matter. 

The number of baryons in a star is the integral over the 
invariant volume element, dei{g^u)d'^x, and the con- 
served baryon number current. The standard textbook 
result is 



No,M = 47r / drr'^no,M{r) 
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2M{r) 



-1/2 



(11) 
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FIG. 2: (Color online) Radii of 1.4 M0 stars (upper panel) 
and maximum masses (lower panel) vs. the relative number 
of mirror baryons to ordinary baryons. 
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FIG. 3: (Color online) Density profiles of neutron stars with 
identical baryon number in the ordinary sector, but different 
number of mirror baryons. The number of ordinary baryons 
is fixed such that M — lA Mq in absence of mirror matter. 



where no {tlm) is the number density of baryons (mirror 
baryons) . The Schwarzschild boundary condition applies 
to the surface at max(i?o, Rm), because the local curva- 
ture of empty space in one sector is affected by the matter 
content in the other sector. The coordinate radius and 
the maximum mass of the sequence depend on the rel- 
ative baryon number in the two sectors, see Fig. [S] In 
Fig. [3] we illustrate the density profiles of NS with identi- 
cal number of ordinary baryons and different number of 
mirror baryons. 

From these results it is clear that the mass-radius re- 
lation of NS is significantly modified in the presence of 
a few percent mirror baryons. The results presented in 
Fig. [2] can be generalized to Nm > Nq by the replace- 
ment M ^ O. It then follows that a mirror NS could be 
accompanied by a compact object in the ordinary sector 
that is only a few kilometers in size. Should they exist, 
such objects would have extraordinary properties, as the 
net gravitational mass is comparable to that of an or- 
dinary NS. The relevance of these results in the mirror 
matter scenario depends on the probability for signifi- 
cant amounts of mirror baryons (baryons) to accumulate 
in the gravitational potential of (mirror) NS and their 
progenitor stars. 

III. CORE COLLAPSE AND LIMITING MASS 

The critical mass of the core of ordinary stars is modi- 
fied in the presence of mirror matter, because mirror par- 
ticles act as an additional source of gravity. Such stars 
therefore collapse at lower core masses and this could re- 
sult in extraordinary compact NS with low mass, due to 
the presence of mirror matter. In order to estimate the 



magnitude of the effect we modify the standard textbook 
calculation of the Chandrasekhar mass of a gravitating 
object supported by electron degeneracy pressure, see, 
e.<?., [13] or [4^. We are interested in cores with high 
density that are near gravitational collapse. The density 
is essentially given by the rest-mass density of nucleons, 
Pq, and the dominant contribution to the pressure, p, 
comes from ultrarelativistic electrons. The EoS is there- 
fore approximately polytropic, p = Kp^, with 7 = |- 
The factor K depends on the relative number of elec- 
trons to nucleons and therefore depends on the chemical 
composition of the core. To leading order, this EoS and 
the hydrostatic equations give a limiting mass of 

M,^1A6Mq(J^^ , (12) 

where = ne/ns is the electron fraction, i.e., the rela- 
tive number density of electrons to baryons. This is the 
traditional expression for the Chandrasekhar mass. 

A self-consistent derivation of the limiting mass re- 
quires that the hydrostatic equations in the two sectors 
are coupled and solved simultaneously, as in Section [ITl 
This is not trivial because the radial dependencies of the 
pressure and density are different in the two sectors. In 
order to obtain the critical mass as a function of the rel- 
ative number of baryons in the two sectors we solve the 
hydrostatic equations, ^ and (pjl- fTU]) . with the poly- 
tropic EoS used in the derivation of the ordinary Chan- 
drasekhar mass. The result is shown in Fig. 21 It is not 
necessary to account for general-relativistic effects in cal- 
culations of the structure of ordinary (and white-dwarf) 
stars, but since the equations are already presented in 
Section |TT] we use them. 
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FIG. 4: (Color online) The modified Chandrasekhar mass 
vs. the relative number of baryons in the mirror sector, 
Nm/{No + Nm), at fixed electron fraction, — 0.5. In 
the limit of 0% mirror baryons the traditional Chandrasekhar 
mass, Mc — 1.46 Mq is reproduced. The critical mass de- 
pends on the electron fraction, i.e., the chemical composition, 
which could be different in the two sectors. For equal electron 
fractions the Ye-dependence of the model is simple, see (|13p . 



We limit the discussion here to the simplified situation 
when the electron fractions in the two sectors are equal. 
In this case the leading-order yg-dependence in the stan- 
dard expression for the Chandrasekhar mass applies also 
to the situation when mirror baryons are present in the 
core, i.e., the numerical solution shown in Fig. |4]is pro- 
portional to (Ye/O.S)^. By fitting a polynomial to the 
numerical results we obtain an expression for the modi- 
fied Chandrasekhar mass 



~ (l.04+1.26g2-l.36g'' + 12.0g6) M©, (13) 

where q = Nm/{No + Nm) — 0.5 is the relative num- 
ber of mirror baryons shifted by an offset of —50% that 
necessarily makes the solution even, because we use the 
same EoS in both sectors. The fit is shown in Fig. |4l 
In the limit of 0% mirror baryons the traditional Chan- 
drasekhar mass, Mc — 1.46 is reproduced. For 50% 
mirror matter our result is consistent with the estimate in 
[i^ l . The nature of the critical mass is somewhat different 
than the picture in this simplified model. In particular, 
the real cause of the collapse is the capture of energetic 
electrons by protons and a softening of the EoS in the 
core. The relativistic mass-limit considered here gives a 
qualitative picture of how the critical mass of the core is 
modified when mirror DM is present, or is accreted from 
the environment. 



IV. ACCRETION OF MIRROR MATTER 

The distribution of mirror DM in galaxies is expected 
to be non-homogeneous, because mirror baryons should 
form complex structures in a similar way as ordinary 
baryons do. There are essentially three different possi- 
bilities for the capture of mirror matter by an ordinary 
NS (the same would apply for the opposite situation, 
i.e., capture of ordinary matter by a mirror NS). First, 
NS should accrete particles from the mirror interstellar 
medium and this must not be in conffict with observa- 
tions, e.f?., by leading to gravitational collapse of young 
NS or too high surface temperatures vs. estimated ages. 
Secondly, the accretion rate could be significantly en- 
hanced if an NS passes through a high-density region of 
space, t.g., a mirror molecular cloud or planetary nebula. 
The third possibility is that NS merge with macroscopic 
bodies in the mirror sector, causing violent events and 
possibly a collapse into a black hole. The probability 
of such events is low and should decrease with the mass 
of the dark objects, but they could be relatively easy to 
discover due to the high energy output, e.g., in the form 
of strange supernova-like events. One such observation 
has recently been reported [1^. We expect that the ac- 
cretion rate increases towards the center of the galaxy 
and where the concentration of dark matter is high, e.g., 
in mirror stellar clusters. Therefore, the effects of mir- 
ror matter should depend on the location and history of 
each star. Note that in general the local environment 
in the hidden sector is different from that in the visible 
sector. For example, an NS located in an empty galactic 
region could be surrounded by stars or a molecular cloud 
in the mirror sector. In fact, the structure formation pro- 
cess is different in the two sectors. Structures at small 
scales, like stars and stellar clusters, are formed essen- 
tially independently in the two sectors because they are 
electromagnetically decoupled. The process of accretion 
of mirror matter in ordinary stars were first studied by 
Blinnikov and Khlopov many years ago [i^, |43| . 

The accretion rate of mirror matter depends on the lo- 
cation of the NS and the structure of the hidden mirror 
sector, which is unknown. We therefore limit the discus- 
sion here to some rough upper-bound estimates, which 
nevertheless leads to a useful conclusion. The accretion 
rate of mirror baryons can be estimated with the result 
for capture of coUisionless particles in the gravitational 
potential of a compact object, Eq. (14.2.21) in 48.] : 

V 10-24 g/cms; \ ^' ^ ' 

Here p^o and Voo are, respectively, the density and speed 
of mirror particles. In the derivation of this expression 
it is assumed that the distribution of mirror particles is 
isotropic and monoenergetic at large distance from the 
NS. The quantities are normalized to the typical values 
for the interstellar medium in the ordinary sector of our 
galaxy. Giant molecular clouds made of ordinary baryons 
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can be tens of parsecs in size and the particle number 
density can reach 10^ cm~'^ in some regions [i^. The 
clouds consist mainly of hydrogen 70%) so the den- 
sity can reach lO"'^^ g/cm'^. If the density of the inter- 
stellar medium in the mirror sector reaches this value, an 
NS at that location would accrete about 10^ kg/s. It is 
possible that mirror molecular clouds have higher den- 
sity than their ordinary counterparts, because the mirror 
sector is subject to different initial conditions and should 
have different chemical composition. In particular, it is 
expected to be helium dominated. We account for these 
uncertainties by increasing the upper-limit estimate for 
the accretion rate from the mirror interstellar medium 
to 10^ kg/s. We have checked that this limit remains 
valid also for NS with arbitrary high kick- velocities. In 
fact, the accretion rate is lower for NS with high speed. 
An accretion rate of 10^ kg/s is, however, far too low 
to have a significant effect on the mass and structure 
of NS. Mirror matter must therefore originate from the 
progenitor star, or from a companion mirror star, if the 
effects on the mass-radius relation discussed in Section 
im are to be observable. Note that the density consid- 
ered above exceeds the average dark-halo density, which 
varies between about a tenth of a GeV/cm'^ to hundreds 
of GeV/cm^ depending on the model used and the loca- 
tion in the galaxy. One GeV/cm^ corresponds to about 
10~24 g/cm'^, so on average the dark halo density is well 
below the 10"^* g/cm^ that we consider in the estimate 
above. 

The gravitational binding energy of NS is of the order 
of 10% of the mass-energy. While the accretion rate is too 
low to significantly affect the mass and structure of NS, 
one could imagine that the associated release of gravita- 
tional binding energy heats the star. The total binding 
energy can be calculated by comparing the energy of an 
NS in equilibrium before and after the accretion of a mir- 
ror baryon. The total mass-energy of a star is M (in units 
where c = 1) and the binding energy is msNB—M, where 
M is the gravitational mass and (Nb) is the baryon 
mass (net baryon number). The change of the binding 
energy when the mirror baryon number is increased by 
one is 
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This quantity is plotted with bold lines in the upper panel 
of Fig. [5] for three different choices of the baryon number 
in the ordinary sector. The EoS used is the same as that 
in Section ini Most of this binding energy is released in 
the mirror sector. The kinetic energy of a mirror baryon 
that falls from great distance to the mirror-matter sur- 
face is ms[l — ^5oo(i?A/)], where gao{RM) is the time- 
component of the metric tensor at the surface. This ki- 
netic energy is plotted with thin lines in the upper panel 
of Fig. [5] and it accounts for more than 97% of the to- 
tal binding energy in all three cases. A significant part 
of the remaining small binding energy is accounted for 
by in- medium binding of the mirror baryon, i.e., tub is 




FIG. 5: (Color online) Change of binding energy per ac- 
creted mirror baryon (upper panel) and the total gravita- 
tional mass (lower panel) vs. the relative number of mirror 
baryons. The baryon number in the ordinary sector is fixed 
such that M{Nm = 0) = 1.2, 1.4 and 1.6 M©. Bold lines in- 
dicate the total change of binding energy psp . which includes 
both gravitational binding energy and nuclear binding energy. 
Thin lines indicate the kinetic energy of a mirror baryon that 
falls from great distance to the surface of the mirror-matter 
core, see text for further information. The solid disks indicate 
maximum-mass configurations, i.e., stars with higher mirror- 
baryon number are gravitationally unstable. 



higher than the energy per baryon in the medium. In 
principle there can be some heating of the ordinary part 
of the NS, e.g., due to weak interactions between the two 
sectors, but as we demonstrate below such effects are ir- 
relevant from an observational point of view. See also 
[sot , where this is discussed in some detail. 

Isolated NS have essentially black-body spectra, so the 
net surface emission is given by the Stefan-Boltzmann 
law, E = AnaT^R^. The observed surface temperature 
of NS is around 10^ K and the radius is of the order 
10 km, so the radiated energy is about 10^^ J/s [52]. If 
we assume that all binding energy is released as heat in 
the ordinary sector (this is clearly an overestimate), the 
upper-limit for the accretion rate, 10^ kg/s, would cor- 
respond to a total heating effect of about 10^^ J/s. The 
thermal emission from the coldest observed NS is about 
two orders of magnitude higher than this upper-limit es- 
timate for the heating. Furthermore, most of the binding 
energy is converted to kinetic energy of mirror particles 
when they fall in the gravitational field. These parti- 
cles will heat the mirror-matter core on impact, but that 
does not affect the temperature of ordinary matter, be- 
cause the two sectors are electromagnetically decoupled. 
Consequently, the accretion of mirror matter in NS does 
not cause significant heating of the ordinary part of the 
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star. The accretion-driven heating should be somewhat 
higher in the case of WIMPs, but the effect is presently 
not observable 

V. CONCLUSIONS 

If DM is made of mirror baryons, they are present in 
all gravitational structures. In particular, they should 
exist in stars at the final stages of stellar evolution. In 
this paper we show that this could have significant effects 
on the formation and structure of NS. Our knowledge of 
NS is limited and highly dependent on observations. The 
situation is complicated by the fact that the conditions 
inside NS cannot be reproduced in laboratories. Ideally, 
we should learn how to calculate the ground-state EoS 
from first principles using the SM and then test it with 
observations of NS, but this task is presently too compli- 
cated. An intermediate goal is to unify "low-density" re- 
sults from nuclear-physics experiments and lattice-QCD 
calculations with observed properties of NS. In that con- 
text it is important to understand the potential effects of 
DM on NS, so that the observations are not misunder- 
stood and false conclusions about the microphysics are 
made. After all, DM apparently is dominant over ordi- 
nary matter and we do not know what it is. In this paper 
we illustrate that the NS sequence is significantly modi- 
fied in the presence of a few percent mirror DM and the 
effect mimics that of other exotic forms of matter, e.g., 
quark matter. An observation of an extraordinary com- 
pact NS is therefore not sufficient to conclude that there 
must be some form of exotic phase of ordinary matter in 
the star, since it can be explained also by the presence of 
mirror matter. 

We find that the NS equilibrium sequence is a func- 
tion of the mirror DM content, i.e., in contrast to the 



equilibrium sequence of ordinary NS it is not a one- 
parameter sequence. The maximum mass of the sequence 
and the radii of stars with typically observed masses, 
M ~ 1.35 Mq, decrease with increasing mirror baryon 
number. The non-uniqueness of the equilibrium sequence 
is the main result of this work, because it is not clear 
that it can be explained with traditional physics, nor 
with self-annihilating DM candidates such as supersym- 
metric particles. New measurements of NS masses and 
mass-radius relations are frequently made, thanks to the 
high resolution and sensitivity of modern observatories. 
The uniqueness property of the NS equilibrium sequence 
could therefore be put to the test in the near future. If it 
turns out that NS masses and radii cannot be explained 
with a one-parameter sequence (after effects of rotation 
are accounted for) , one could search for a correlation be- 
tween the orbits of the observed NS and the expected 
density of DM along those orbits. The accumulation of 
mirror matter in stars could have a significant effect also 
on the critical mass for core collapse. This is a poten- 
tial formation mechanism for NS with low masses and 
extraordinary small size. These results can be used in 
future studies to test models of the mirror Universe. 
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